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The existence of large-scale magnetic fields in galaxies is well established, but 
there is no accepted mechanism for generating a primordial field which could grow 
into what is observed today. We discuss a model which attempts to account for the 
necessary primordial field by invoking a pseudo-Goldstone boson coupled to electro- 
magnetism. The evolution of this boson during inflation generates a magnetic field; 
however, it seems difficult on rather general grounds to obtain fields of sufficient 
strength on astrophysically interesting scales. 
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I. INTRODUCTION 



The existence of a magnetic field of ~ 10 -6 gauss in our galaxy and other galaxies is well 
established 0. The explanation of how such a magnetic field arose is, however, far from 
certain. While the creation and evolution of stellar magnetic fields is fairly well understood, 
the extension of these theories to galaxies suffers from problems relating to both scale length 
and time scales. 

Zeldovich, Ruzmaikin, and Sokoloff M and Parker |2| discuss the origin and effects of 
magnetic fields in the universe and draw the conclusion that the galactic field arises from 
a dynamo mechanism. The dynamo model requires a seed field at the epoch of galaxy 
formation which is coherent over a scale of ~ 1 Mpc. We can parameterize the strength 
of a primordial field by r = ps/p 7 , the ratio of the energy density p B = B 2 /8n in the 
magnetic field relative to that of the background radiation p 7 . (This ratio is constant while 
the universe is a good conductor, which is almost always 0.) Then the field required to 
seed a galactic dynamo satisfies r > 1CT 31 , corresponding to an intergalactic field at the 
epoch of galaxy formation (z ~ 3-5) of <^ 10~ 20 gauss. The implications of this requirement 
in terms of the origin of such fields and their possible effect on star formation (or the early 
history of the galaxy in general) are discussed by Rees [Q. 

Kulsrud has argued that the galactic dynamo explanation is fundamentally flawed 
in that the mean square deviations of the magnetic field will grow much faster than the 
mean field itself, resulting in a disordered field with a much smaller mean field strength 
than would naively be expected. Kulsrud then discusses the possibility that the magnetic 
field originated in the early universe and was embedded in the medium out of which the 
galaxies ultimately formed. In this case, Kulsrud argues that the intergalactic field would 
need to be ^ 10 -12 gauss at the epoch of galaxy formation (giving r > 0~ 15 ) in order to 
account for the observed interstellar field. 

Any dynamo theory requires a mechanism for generating the required seed field; how- 
ever, compelling mechanisms have been elusive. The hot plasma in the early universe is 
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highly conducting and thus should strongly inhibit the growth of a magnetic field, even to 
r ~ 10~ 31 . Furthermore, any hypothetical process at work in the very early universe must 
be able to produce fields with characteristic length scales much larger then the horizon at 
that time, in order to correspond to galactic scales today 

The advent of inflation |Bj has opened the door to new possibilities for generating a 
primordial magnetic field. There are two key features of an inflationary universe that make 
the possibility of creating a magnetic field during this time particularly attractive. First, 
if there were an inflationary epoch at very early times in the universe, the exponential 
expansion would have reduced the conductivity to a negligible value by reducing the charged 
particle density, thus allowing the creation of a substantial magnetic field. If this field were 
then frozen into the plasma created during the subsequent reheating of the universe, it will 
be supported by the effectively infinite conductivity of the plasma so that its strength will 
decay only as inverse square of the scale factor. 

Second, if inflation did occur, then the entire observable universe today was, at some 
point in the early universe, contained entirely within the particle horizon. It would then 
be possible to use physical mechanisms operating on a scale smaller than the horizon to 
generate magnetic fields that are coherent over macroscopic scales today, an opportunity 
which is not available in models of the early universe without inflation. 

Turner and Widrow || (TW) investigated the possibility that quantum fluctuations 
during an inflationary epoch might have generated a magnetic field that could be sustained 
after the wavelength of interest crossed beyond the horizon and thus give the observed field 
today. TW considered coupling the electromagnetic field to the curvature tensor so as to 
amplify the fluctuation-induced field, but found satisfactory results could be obtained only 
at the expense of breaking gauge invariance. 

Ratra has argued that it is possible to generate a magnetic field with present field 
strength of ^ 1CT 10 gauss on a scale of 1/1000 the present Hubble radius by coupling a 
scalar field <3> to the electromagnetic potential through a term of the form e^F^F^ , 
where F^ = d^A u — d u A^ is the electromagnetic field strength. This would be sufficient to 
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explain the galactic magnetic field, but it remains to be seen whether or not this coupling 
could arise naturally in realistic particle physics models. 

TW also suggested that the magnetic field could be sustained by coupling the EM field 
strength to a pseudoscalar axion field <fi via an interaction term in the Lagrange density of 
the form 



where F^ LU = \t^ vtxJ F pa is the dual of F^, and ^ 77 = (a/2n)/f, where / is a coupling 
constant with units of mass and a is the fine structure constant. However, they did not 
complete the necessary analysis to show whether or not this might indeed be the case. 

An interaction of this form has been studied by Carroll and Field || (see also ||), who 
found that the evolution of a Fourier mode of the magnetic field with wavenumber k is 
governed by 



where F± = a 2 (B y ± iB z ) (the ± refers to different circular polarization modes of the 
magnetic field), di] = dt/a, and a is the scale factor of the universe (normalized so that 
ao — 1 where ao is the value of a today). One (or both) of the polarization modes will be 
unstable for k < g^^dcf)/ dr]\, where both polarization modes can be unstable to exponential 
growth if is oscillatory. Thus, if such a scalar field exists during inflation (perhaps the 
inflaton itself (1(||) with the above coupling, this might provide a mechanism for generating 
a substantial magnetic field. 

Here we will consider a generalization of the possibility suggested by TW ||, coupling 
the photon to an arbitrary pseudo-Goldstone boson (PGB) rather than the axion of QCD. 
The PGB is characterized by a spontaneous symmetry breaking scale / (as above) and a 
soft explicit symmetry breaking scale A (see Sec. II). We find that significant growth occurs 
only at a temperature near A, and that the magnetic field strength thus generated cannot 
give an astrophysically interesting field at the end of inflation. 



C = ^F^ , 



(1.1) 




(1.2) 
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We should point out that, while our notation throughout this paper suggests that we are 
working with the photon of the standard U(l) em symmetry, the photon as a separate U(l) 
gauge boson will not exist at high temperature, since the SU(2) x U(l) gauge symmetry 
will not have been broken. Nevertheless, our results should be correct up to factors of 
order unity simply because the U(l) hypercharge symmetry projects onto the photon with 
a multiplicative factor of cos9 w 0.88 at the electroweak phase transition. 

We will use units in which h — c — ks — 1, such that G = m~ 2 , where m p i « 1.22 x 10 19 
GeV is the Planck mass. 

II. THE SET-UP 

In this section we briefly review the essentials of inflation, as well as the physics of 
pseudo-Goldstone bosons and their couplings. 

Throughout this paper, we will assume that the universe is in a spatially flat FRW 
cosmology in which the metric is given by 

ds 2 = a 2 (r]) {-drf + rfx 2 ) , (2.1) 

where x represents the standard Cartesian 3-space (comoving) coordinates. In addition, we 
will assume that the universe is a perfect fluid with the equation of state p = jp, where 
p is the pressure, p is the total energy density, and 7 is a constant. Using this equation 
and energy-momentum conservation, it is straightforward to show that p oc a~ 3 ( 1+ ^ which, 
from Einstein's equation, gives a oc t 2 / 3 ( 1+ T). In order to explain the horizon problem, we 
require that, at some time in the past, the scale factor was growing faster than the horizon 
(H -1 , where H = a/a is the Hubble parameter and an overdot denotes differentiation with 
respect to physical time t). Thus, we require — 1 < 7 < — |. For simplicity in the following 
discussion we will restrict ourselves to inflation in which 7 = — 1. This value for 7 gives the 
best possible conditions for generating a magnetic field, since the amount of inflation from 
the time that a given comoving wavelength crosses outside the horizon is minimal in this 
case, so this does not limit the validity of our results. 
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At the end of inflation, the universe enters a reheating phase, in which the energy density 
is matter-dominated. As a simplifying assumption, we take the process of reheating to be 
instantaneous, such that the universe goes directly from inflation to radiation domination. 
Once again this is a best-possible assumption, since the magnetic field will decay more 
rapidly (relative to the total energy density) during a matter-dominated phase, in which 
p oc a~ 3 . 

Standard inflation is characterized by two parameters: the mass scale for the total energy 
density M = p 1//4 (note that this is a constant since p is constant during inflation with 
7 = —1); and the temperature Trh to which the universe reheats at the end of inflation. H 
is then given by 

Air 

H 2 = -^Vm 4 . (2.2) 

If we assume that the universe expands adiabatically after inflation so that the entropy per 
comoving volume element remains constant, it can be shown (see e.g. ||11|| ) that the total 
expansion from the time a given comoving wavelength A (which is equal to the physical 
wavelength today due to the normalization of a) crosses outside the horizon until the end 
of inflation is given by 

i 



a inf _ 1n 26 X [M 2 TrhY 

l pl 



- 10 26 -- -^A^- , (2.3) 
d\ Mpc I ml ' 



where ai n f is the value of a at the end of inflation (but before reheating), and a\ is the value 
of a at the time A crosses outside the horizon. Expressing this in terms of the number of 
e-foldings N\ (a in f/a\ = e Nx ) in the expansion, we have 

iV A ~48 + lnf-M+glnf f U-lnf ? R h \ ( } 

\MpcJ 3 UO 14 GeW 3 UO 14 GeV J v ; 

Our assumption that reheating lasts for a negligible time amounts to setting M = Trh. 

In this paper we are concerned with the pseudo-Goldstone boson of a spontaneously 
broken symmetry. PGBs are characterized by two mass scales: a large mass / at which the 
global symmetry from which the PGBs arise is spontaneously broken, and a smaller scale 



A at which the symmetry is explicitly broken. For concreteness, we imagine the breakdown 
of a global U(l) symmetry, resulting in the familiar Mexican hat — the radial degree of 
freedom gets a vacuum expectation value of order /, and the angular degree of freedom 
becomes a massless boson <fi. The hat is tilted by a small term of order A; the formerly 
massless scalar becomes a PGB with mass of order 

m^A 2 /f. (2.5) 



Cosmological constraints on the parameters / and A have been studied in |L2 . 

In many models, PGB's interact with fermions by coupling to the axial vector current 
(for a review of PGB's and their couplings, see fl3f): 

Ant = -Ut"75^ = AW ; (2-6) 

well-known examples include pions and axions. Since the symmetry associated with the axial 
current (that is, chiral symmetry) is anomalously broken, the current itself is not conserved: 

= ^-F, U F^ , (2.7) 
Ztx 

where a is the fine structure constant. Integrating by parts, we find that the anomaly (2.7) 
induces a coupling between and the electromagnetic field of the form (1.1). 

A similar situation occurs in the low-energy limit of string theory, which involves an 



antisymmetric two-index tensor field B pL , [13,14]. The Lagrangian for Bn U includes a kinetic 



term H^ l/p H' J ' up , where H^p is an antisymmetric field strength tensor. The demand that 
the theory be free of anomalies requires that the definition of H pvp include a term involving 
gauge bosons (which we take to be abelian for simplicity): 

Hp,v P = d[p,Bvp\ ~ \pPvp\ > (2-8) 

where square brackets denote antisymmetrization. In four dimensions the equations of mo- 
tion for B pv allow us to recast its dynamics (at least semiclassically) in terms of a pseu- 
doscalar <fi by the identification 



a/3-f 



(2.9) 



The interaction between this pseudoscalar and electromagnetism, as implied by (2.8) and 
(2.9), can be described by an effective Lagrangian of the form (1.1). Thus, string theory 
offers the possibility of PGB's of the type we discuss. 



We envision a scenario in which the desired magnetic field is entirely created during 
inflation and then frozen into the post-reheat plasma, so that the field strength decays as 
the inverse square of the scale factor after reheating. Thus, under the assumption that the 
dominant component of background photons is created at reheating, both ps and p 7 will 
decay as a -4 , such that r = psj 'p 7 is essentially constant after inflation. As mentioned in 
the introduction, we seek r = r > 10 -31 (to seed a galactic dynamo) or > 1CT 15 (to directly 
account for the observed magnetic field) at scales of ~ 1 Mpc today. 

For convenience we define r during inflation as the ratio of pb to the total (vacuum- 
dominated) energy density. From above, we can readily calculate the required value of r at 
the time that the wavelength of interest (comoving scale 1 Mpc) crosses outside the horizon 
in order to generate astrophysically interesting fields today. This is given by 



Since r < 1, we can use this to calculate the maximum allowed value for M such that 
there would be any hope for meeting the above condition. Using the values given above 
for r , this gives M < 10 GeV for a galactic dynamo, or < 1 MeV to directly account for 
the field. It is important to realize that these are extreme upper bounds that can only be 
obtained under ideal situations (i.e. the magnetic field energy being comparable to that in 



III. BASIC REQUIREMENTS 




(3.1) 



or, using (2.3), 




(3.2) 
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the vacuum at the time the wavelength of interest crosses the horizon, and with negligible 
reheating). Even so, the upper bound for seeding the galactic dynamo is, at best, marginal 
(that is, while current constraints on the time at which inflation might occur could allow it 
to occur as late as a temperature of ~ 1 GeV since we merely require the universe to be 
radiation dominated by the epoch of primordial nucleosynthesis ||, constraints arising from 
baryogenesis will probably require M > 200 GeV, at least). Thus, even if the magnetic field 
has an energy density comparable to the background energy density when 1 Mpc crosses 
the horizon during inflation, it will be too weak to explain the observed magnetic field if its 
energy decays according to the normal a -4 . 

In other words, in order to generate a significant magnetic field during inflation, we 
require "superadiabatic growth" — that is, a mechanism that will continue to increase the 
energy density (or, at least, decrease the rate of decay) of the magnetic field at a wavelength 
of 1 Mpc after 1 Mpc becomes super-horizon sized. At first sight, this may seem impossible 
due to the fact that such a mechanism must apparently act ill cLll ctCctUScll way. However, it 
is possible that inflation may create a field that is coherent over scales much larger than the 
horizon, and that this field can subsequently generate a magnetic field that is also coherent 
at super-horizon scales simply by classical field interactions. 



IV. EQUATIONS OF MOTION 

The Lagrange density for the photon and scalar field is 

C = - Qv^W + V{(j>) + \f, u F^ + ^F^F^ , (4.1) 

where g = —det^g^) and V M denotes the covariant derivative. As mentioned previously, we 
are considering only the U(l) fields and ignoring any possible effects from the nonabelian 
gauge fields. 

The equations of motion for are 

- V^V,0 + ^ = - 9 -fL F , u F^ , (4.2) 
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and the equations of motion for F^ are 

= - 9<hl (V^)F^ , 

along with the Bianchi identity 



V M F^ = . 



The equations of motion are more transparent if we define the E and B fields by 



(4.3) 



(4.4) 



= a 



Then (4.2) becomes, after expanding the covariant derivatives, 
d 2 (J) d4> _ , odV((f)) 



E x 


Ey 


E z \ 





B z 


-By 


-B z 





B x 


By 


-B x 


o ) 



(4.5) 



+ 2aH^- - V 2 + a 2 — ^p- = g^a 2 E ■ B , 



(4.6) 



drf dr\ d(f> 

where V represents the usual 3-space gradient (for comoving coordinates). Similarly, (4.3) 
becomes 

d 



with 



^-(a 2 E) - V x (a 2 B) = -g Hj ^ a 2 B - ^ 77 (V0) x a 2 E , 



V • E = - g<h7 (V<j>) ■ B , 



(4.7) 



(4.8) 



while the Bianchi identity becomes 

d 



a 2 B) + V x (a 2 E) = , 



(4.9) 



with 



V-B = 



(4.10) 



Since we are interested in the specific case where the background space-time is inflating, 
we make the assumption that the spatial derivatives of are negligible compared to the 
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other terms (if this is not the case at the beginning of inflation, any spatial inhomogeneities 
will quickly be inflated away and this assumption will quickly become very accurate). Then, 
eliminating E in the above equations, we have 

(^-V 2 -^|vx)(a^B) = 0. (4.11) 

Taking the spatial fourier transform of this equation so that 

B(r ] ,k) = ^j e* k - x B(7?,x)d 3 x, (4.12) 
and writing F = a 2 B, we then have 

+ A; 2 F-^ 77 ^kxF = O. (4.13) 
Finally, if we take k to point along the x-axis and define F± = F y ± iF z , this becomes 

^ + ^± 9tll d A k y ± = . (4.14) 

We can similarly manipulate the equations in an attempt to produce an expression for the 
evolution of E. It turns out, however, that we cannot uncouple the E field from the B field. 
We have, in short, 

(J^ - V 2 - ff^Vxj (a 2 E) = -^ 77 0a 2 B , (4.15) 

which, after taking the space fourier transform and defining G = a 2 E and G± = G y ± iG z , 
becomes 

d 2 G ± f l2 d(j)\„ d 2 <f) „ , 

+ \k ± 9<hr^jk ) G ± = -9^^F± . (4.16) 

In order to determine the evolution of E and B, we need to know how evolves. We look 
first at the case where E and B make a negligible contribution to the equation of motion for 
(j). Furthermore, we will consider the evolution after the time when the explicit symmetry 
breaking for the PGB becomes important (at a temperature scale A). The potential for the 
angular degree of freedom in a tilted Mexican hat is 
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\/(0)=A 4 [l-cos(0//)] . (4.17) 

Since the details of the potential do not affect our results, we will expand to lowest order: 
V(4>) ~ (A 4 /2/ 2 )0 2 . We again ignore spatial derivatives in 0, as well as any back reaction 
from the electromagnetic fields, to give 

where we have used dr] — dtj a to write this in terms of physical time rather than conformal 
time. The general solution to (4.18) will be approximately 



<f)(t) « / exp 



~\ (3#±yW-4^)(*-* ) 



oc < 



a -3/2 sin /Al t for Ai >jf J ) 

4 ; a ( 4 - 19 ) 

exp y-^pt) for y < 



where we have used the fact that a oc exp Ht. 

Looking again at (4.14) we see that we will only have a growing mode for the magnetic 
field if 

~T~ = a ~T > k ■ ( 4 -20) 
drj dt v 1 

Furthermore, since we want this growth to be as large as possible, we will choose k such that 
k ~ g^a&max at the time of interest. For our initial analysis, we will assume A 2 // ^> H so 
that is oscillating rapidly compared to changes in a. Also, this implies that a is essentially 
constant over several oscillations in 0, and thus we can write Arj ~ At /a where a is constant 
for time intervals At ~ //A 2 . 

In order to estimate the total growth in F±, we note that, for a fraction e (where e is not 
necessarily small, although numerical integration of these equations for some cases indicates 
that e ~ 0.1) of each period, we can write F± oc e aAri where a = yj g^kdcf) / drj ~ ag^A 2 and 
Ar] At /a ~ ef /aA 2 . Furthermore, this will continue for a time ~ H^ 1 (since this is the 
time scale on which a and the amplitude of are changing), or for a total of n ~ H~ x A 2 jf 
oscillations, from which we can estimate the total growth in F± as 
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|k „ exp . (4.21) 

Of course, the exponent may contain other factors of order unity, but this estimate allows 
us to understand the dependence of the growth in the magnetic field on the parameters in 
the problem. Note that the analysis is similar if A 2 // ~ H, but then Arj w At /a where a 
is constant is only true for time intervals At < H~ l . So using Arj ~ if -1 /a, we see that we 
get exponential growth with the same expression in the exponent, but now this exponent is 
much smaller. 

Since we are interested specifically in long wavelength magnetic fields, we would also like 
to know the wavelength at which we get the most growth by this mechanism. This follows 
from our assumption that the maximum growth occurs at k ~ ag^^max ~ ag0 77 A 2 , which 
gives 

A » , (4.22) 

where A is the wavelength today. More importantly, at the time the growth occurs the ratio 
of the wavelength to the horizon length is given by 

aX 2naH 2nH 



H 1 k S^ 77 A 2 

But, looking back at (4.21), we see that this ratio is, essentially, just the inverse of the factor 
in the exponent, i.e. the larger the growth in the magnetic field, the smaller the wavelength 
at which it occurs! Further, since an increase in the wavelength by a factor (5 will result in 
a decrease in the exponent for the growth by a factor f3~ 1 ^ 2 , it is apparent that significant 
growth in the magnetic field will only occur for wavelengths that are sub-horizon sized. 

We have shown that we cannot have growth in the magnetic field for super-horizon sized 
wavelengths when the back effects of the E and B fields on are small, but there is still 
the possibility that the interaction between the fields and could allow the magnetic field 
to be sustained such that it does not decrease as a~ 2 . However, from (4.7) and (4.8), we 
see that, if decays as a~ 3//2 (that is, behaves as a non-relativistic fluid), then the RHS's 
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of these equations will rapidly become negligible compared to the individual terms on the 
LHS's. We then recover the source-free Maxwell equations, from which B still decreases as 
a~ 2 . (The situation is only exacerbated if behaves as a relativistic fluid.) In short, there 
seems to be no way of sustaining the magnetic field at super-horizon sized wavelengths. 

V. SUMMARY 

In this paper, we have considered a mechanism for creating a large scale magnetic field 
during inflation, proposed by TW, in which the magnetic field is coupled to a pseudo- 
Goldstone boson. We showed first that the scale required of the magnetic field in order to 
explain the galactic magnetic field (~ 1 Mpc) implied that the growth had to occur at super- 
horizon sized wavelengths since the uncoupled equations of motion for a U(l) gauge field 
imply that the energy density in the field would simply decay too quickly to be significant 
at the end of inflation if there was no enhancement in the field for wavelengths larger than 
the horizon. 

We then considered the classical evolution of a U(l) gauge field coupled to a PGB when 
the back-effects of the field on the PGB were negligible and showed that such a coupling 
can, in fact, produce growth in the field. However, such growth can occur only at sub- 
horizon wavelengths, and thus does not provide a solution to the above problem. In the 
more general case when we allow the back-effects of the U(l) field to be significant, we still 
have the problem that any natural decay of the field ultimately allows the U(l) field to 
uncouple from the PGB leaving us once again with a free U(l) field. Hence, it seems to be 
impossible to create a significant magnetic field by simply coupling the magnetic field to a 
PGB during inflation. 
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